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We study the properties of the B — > tt and B — > K transition form factors in partially quenched 
QCD by using the approach of partially quenched chiral perturbation theory combined with the 
static heavy quark limit. We show that the form factors change almost linearly when varying the 
value of the sea quark mass, whereas the dependence on the valence quark mass contains both the 
standard and chirally divergent (quenched) logarithms. A simple strategy for the chiral extrapola- 
tions in the lattice studies with N se& = 2 is suggested. It consists of the linear extrapolations from 
the realistically accessible quark masses, first in the sea and then in the valence quark mass. From 
the present approach, we estimate the uncertainty induced by such extrapolations to be within 5%. 
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INTRODUCTION 



In order to extract the Cabibbo-Kobayashi-Maskawa 
matrix element \V u b\ from the experimentally measured 
decay rate for B —> ixtvg at Belle and BaBar, a reli- 
able theoretical prediction of the corresponding hadronic 
matrix element is indispensable. Two major sources of 
systematic uncertainty in the current lattice QCD calcu- 
lations of the heavy — > light pseudoscalar meson tran- 
sitions are the quenched approximation and the errors 
associated with the chiral extrapolations. All the avail- 
able results of the lattice studies of the semileptonic 
B — > 7r decay are obtained in the quenched approxima- 
tion (N = 0) and by working with light quark masses 
larger than m p s hys /2 

Since fully unquenched QCD lattice simulations (i.e., 
with N = 3) are not feasible it is important to have a 
method to assess whether or not a given physical quan- 
tity is prone to large quenching errors. With such an 
ambition in mind, Sharpe pj, and later Bernard and 
Golterman p|, formulated the quenched chiral pertur- 
bation theory (QChPT). By confronting the predictions 
derived in QChPT with those obtained in the full chi- 
ral perturbation theory (ChPT), one gets a rough es- 
timate on the size of quenching errors. The approach 
has been extended to the heavy-light quark systems by 
combining ChPT and the heavy quark effective theory 
(HQET) 0,0. The effect of complete quenching on the 
lattice determination of the B — > tt form factors has been 
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studied recently by the present authors in Ref. 7. It has 
been found that the quenching errors on the B — > tt and 
B — » K transition form factors may be uncomfortably 
large (typically larger than 20%). This conclusion some- 
what spoils the significance of the impressive agreement 
that has been reached amongst various lattice groups us- 
ing different lattice techniques to compute B — > tt decay 
in the quenched approximation [l[. 

It is therefore highly important to perform the simu- 
lations in which the effects of dynamical quarks are in- 
cluded. A step in that direction is to implement the 
partial quenching, i.e., to include in the simulation 7Vs ea 
dynamical light quarks with masses generally different 
from those of the valence quarks. The lattice studies 
with iV sea = 2 are likely to be performed first, which is 
the main motivation for the present work. 

We calculate the form factors for B — > tt and B — > K 
transitions in partiall y qu enched chiral perturbation the- 
ory (PQChPT) HllEa, combined with HQET 6]. We 
work in the static heavy quark limit and at next-to- 
leading order (NLO) of the chiral expansion. This ap- 
proach is valid for small recoil momenta (v-p), i.e., the 
same one currently accessible from the lattice simula- 
tions. We limit ourselves to the case of N sea degenerate 
sea quarks of mass m sea . On the basis of our calculation 
we conclude that the dependence of the form factors on 
the sea quark mass is essentially linear, and that the form 
factors are finite as m sea — >• when m va i. is non-zero. On 
the other hand, the limit m va i. — > 0, with m sea ^ ra va i., is 
not well defined, since in this case the form factors con- 
tain the chirally divergent "quenched" logarithmic terms 
m sea lnm va i.. Our analysis of PQChPT with N sea = 2 
shows that a simple linear chiral extrapolation first in 
the sea quark mass, m sea , and then in the valence quark 
mass, m va io introduces an extrapolation error of only 



2 



5%, where the linear extrapolations are made from the 
range of the quark masses that are currently accessible in 
the lattice simulations. Furthermore, by assuming that 
the low energy constants in the full ChPT with N = 2 
and N — 3 are equal, we deduce the residual quench- 
ing errors of lattice studies with 7V sea = 2 to be in the 
ball park of 10 — 20%. Those two conclusions indicate 
that, in comparison with the fully quenched simulation, 
the additional computational cost of making the partial 
unquenching (with N sea — 2) is well worth the effort. 

The remainder of this paper is organized as follows. In 
Sec. |n]we remind the reader about the elements of the 
PQChPT and set our notation. In Sec. lIIII we present the 
results of the NLO calculation in PQChPT for B -> tt(K) 
form factors. The quenching errors and chiral extrapola- 
tion formulas are discussed in Sec. IIVI Main conclusions 
are shortly stated in Sec. 



II. PARTIALLY QUENCHED CHIRAL 
PERTURBATION THEORY 



In this section we make a succinct summary of 
PQChPT. We consider partially quenched QCD with 
7V sea sea quarks degenerate in mass (m sea ), and with two 
valence quarks q a and q b of masses m a and ra&, respec- 
tively. Partial quenching is introduced by adding bosonic 
"ghost" -quarks q a ^ of spin \ and mass which cancel 
the fermion loops of valence quarks @, Lj U ] Assum- 
ing the spontaneous symmetry breaking pattern to be the 
same as in full QCD [l£j, the leading-order Lagrangian 
for the (pseudo) Goldstone bosons is [3, 0, EU 



ght 



^ str (d^d^) + ^ str(M£ + 
8 2 

+A, (1) 



with / « 130 MeV, str (17) = ^Uue h and 



M 



diag(l,..., 1,-1,-1), 

iVsea + 2 

diag(ra a , m bl m sea , • • • , m sea , m a , m b ), 



while 



exp ( 2ij 



x} 
x 4> 



(2) 



(3) 



The graded SU(2 + 7V sea |2) matrix $ contains a [(2 + 
N se &) x (2 + iV sea )] matrix <j) of quark-anti-quark bosons, 
a [2 x 2] matrix <j) of ghost-anti-ghost bosons, and the 
matrices of pseudoscalar fermions, ~ QQ an d X ~ QQ- 
In this framework the heavy rj' has been integrated out 
(mo — > oo ) in a way similar to the full ChPT. This is to 
be contrasted to the fully quenched theory, where the rj f 
cannot be integrated out. The propagators of the "off- 
diagonal" and "diagonal" mesons (^> a 6 ~ q a Qb) are given 



respectively by 

($ab\$ba)a^b 
($aa\$bb) =i 



Sab^a 1 



- M* a N sea (p* - M 2 J (p* - Ml h) 



(4) 



with M 2 6 



2jj j0 (m a + m h ), and Mj 4fi m se3i . 
Since PQChPT with 7V sea = 3, and m a = 7715 = m sea 
coincides with the full ChPT with m u = rrid = ra s , the 
corresponding low energy constants obtained in the two 
theories are equal. 1 Since the low energy constants do 
not depend on the quark masses, the mentioned equality 
among low energy constants persists even when the va- 
lence and sea quark masses of partially quenched theory 
are not of the same size On the other hand the low 
energy constants do depend on the number of dynamical 
quark flavors and therefore the equality holds only for 

Nsea = 3. 

The two terms appearing at NLO in the effective La- 
grangian JU relevant to B — > tt and B — > K form factors 
are 

U = 4/i {^4str (^S^S^str (A^Et + Y>M*) 

+L 5 str [d^d^Z (M^ + ZM^)] +•••} .(5) 

It is straightforward to verify that Eq. leads to the 
standard, full QCD, chiral Lagrangian after setting str — » 
tr, <1> — » (j) (for a review of ChPT see Ref.[l2|or any paper 
listed in Ref. E|). 



A. Incorporating the heavy quarks 

To extend the applicability of the Lagrangian JJJ to 
the heavy-light mesons, the heavy quark spin symmetry 
needs to be included. This is achieved by combining the 
pseudoscalar (P a ) and vector (P* a ) heavy-light mesons 
in one field: 



H a (v) 



H a (v) = 7oi?I(w)7o • 



(6) 



We also introduce the covariant derivative and the axial 
field as 

D b ;H b = d^H a -H b h^d^ + ^] ba , 



(7) 



1 Note, however, that in the usual definition of low energy con- 
stants 11] an SU(3) relation between operators is used IT3 . 
which is not valid for the general SU(A^i|A^2) case. This subtlety 
does not concern us here, as it does not change the definition of 
low energy constants entering in our calculations. 
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where a and b run over the light quark flavors, and £ = 
a/S. Finally, the partially quenched chiral Lagrangian 
for the heavy-light mesons in the static heavy quark limit 
reads 0,0 

£heavy = - str a Tr[H a ivD ba H b } 

+g str a Tr [H a H b ^ A^ 75 ] + £3 , (8) 

where g is the coupling of the heavy meson doublet to 
the Goldstone boson. The higher order terms in the 
expansion in v-p and in m q [0(p 2 )] have the following 
form HQ: 

£3 = 2\istr a Tr[H a H b ](M+) ba 

+fci str a Tr[ff it;.2) teJ ff6](^+)ca 

+fc 2 str a TT[H a ivD ba H b ] str c (M+) cc + . . . ,(9) 

with M+ = 7}(^M^ + £M£). We again display only the 
terms that contribute to the heavy-to-light form factors 
of B — > 7r, if transitions. In the above equations, "Tr" 
stands for the trace over Dirac indices, whereas "str" is 
the super- trace over the light flavor indices. The stan- 
dard chiral Lagrangian for heavy-light mesons 0, 13 * s 
recovered by replacing str — > tr, <I> — > 0. 

The bosonized fteaw/ — > weak current (V — A), 
in the static heavy quark limit and at NLO in the chiral 
expansion, reads [5|, Q 



i) 



1 C\ 1 

J" = g o7 "(l - 75 )Q -> — Tr[ 7 ^(l - 75)^4 



T 



X! Tr[ 7 ^(l - 7s)ff c ]4(^+: 



c6 



(10) 



x 2 Tr[ 7 ^(l - 75)^]4 str c (M + ) cc . 



The phase of the heavy meson can be chosen in such 
a way that the constants a, K\ and x 2 are real. At the 
leading order in the chiral expansion, the constant a is 
proportional to the heavy-light meson decay constant, 
a = ^rriufB- 



B. Form factors 



B 



B 



FIG. 1: The point 1) and the pole 2) tree level Feynman di- 
agrams contributing to heavy — > light transition form factors. 
The box denotes the weak current insertion. 



We work in the static limit (ra# —> 00) in which the 
eigenstates of QCD and HQET Lagrangians are related 
through —L=\B(p B )) qcd = |#0))hqet. In the static 
limit it is more convenient to use a definition in which the 
form factors are independent of the heavy meson mass. 
We will use 



(P(p)\q^b v \B(v)) mET 



[PtJL ~ (v-p)v„]f p (v-p) 
+ vMv-p) , (12) 



where the field b v does not depend on the heavy quark 
mass. The form factors f PiV are functions of the variable 



v-p 



2m b 



(13) 



which is the energy of the light meson Ep, in the heavy 
meson rest frame. The relation between the quantities 
defined in Eqs. ([TT]) and is obtained by matching 
QCD to HQET at the scale \i ~ m b [l6|- By setting the 
matching constants to their tree level values and neglect- 
ing the subleading terms in the heavy quark expansion, 
one has 



Fo(q 2 ) 



-fv(v-p) 



-f P (v-p) ■ 



(14) 



i.e., the usual heavy mass scaling laws for the semilep- 
tonic form factors [l7|. 2 



The standard decomposition of the vector current ma- 
trix element between two pseudoscalar meson states is 



(P(p)\qi»b\B(p B )) 



f \ 1 

VPB +P)n - q^- 



+- 



rn. 



rat 



III. EXPRESSIONS FOR THE FORM FACTORS 
IN PQCHPT 



In this section we give the expressions for the form 
factors f p (v-p) and f v (v-p) as derived in PQChPT with 

(11) 



where the form factors F +5 o are functions of the momen- 
tum transfer squared q 2 = (pb — p) 2 • A light meson P 
stands for 7r, K with the light quark in the current being 
q = d,s, respectively. 



2 Notice that in the heavy quark limit the tensor current 
form factor, Ft((? 2 ), defined as (P(p)\qa^ Liy q u b\B(pB)) = 
i [q 2 (PB +P)v - {m 2 B - m 2 p )q^] F T (q 2 )/(m B Imp), is related 
to the vector form factor via the Isgur-Wise relation lilt . 
F T (q 2 ) = F+(q 2 ). 




A/sea degenerate quarks of mass m sea , and with the va- 
lence light quarks of masses m a and m&. When necessary, 
we will use to denote the light pseudoscalar meson 
with the valence quark content P^ a ~ q^q a and B a ~ bq a 
for the heavy mesons. 

The tree level expressions for B a —> P^ a transition form 
factors are given by the point and pole diagrams in Fig.d 
which give rise to f v and f pi respectively 



where A J = ra#* — ^B a - Although the heavy quark spin 
symmetry suggests — > 0, we keep A^ finite in the 
tree-level form factors because it provides the pole to the 
form factor f p at m B * . 3 

The NLO chiral corrections to the form factors are 
conveniently expressed as 

f P B r Pba = fZ ee { 1 + ( s fp,v) Loop 

+ c^m a + 4' v m b + & sea ], (16) 



3 The pole dominance is easily seen if one rewrites the denomina- 
tor of f P (v-p) as v-p + A* = (tub* /2) (l - g 2 / m B*)' where the 
corrections (m^* — m^)/m^ and rap/ra^ are neglected. 



where 

<=cl-*i = \k 1 - 1 -^, (17a) 

C 6 = C b + ^1 = x l ' ( 17b ) 

<4a = C =(^2 + \k 2 - 3 ~^y sea , (17C) 

that arise from the 0(p 3 ) (0) and (9(p 4 ) terms in the 
Lagrangian, as well as from the weak current (|TU|) . The 
loop contributions are written as 

(Sf P , v ) Loop = 5 fp'v + \ Sz b OOP + \ sz p° OP - (18) 

where the sum runs over all the graphs depicted in Fig.0 
and the last two terms arise from the loop contributions 
to the wave- function renormalizations. The explicit ex- 
pressions for the loop corrections in Eq. (|18|) are rather 
lengthy, and we relegate them to Appendix |£l In the 
calculation of the loop integrals we used the naive di- 
mensional regularization, and the MS + 1 renormalization 
prescription, i.e. we subtract A = 2/(4 — d) — 7 + ln(47r) + 
1 [Tl|. We neglect mass differences A between P a , P5, 
P* and Bl meson states whenever they appear in loops. 
We also remark that no dependence of the form factors 
f v and f p on v-p arises from the counterterms, so the 
modification of the tree level v-p dependence is entirely 
due to the chiral loop corrections. 

Finally, we explicitly checked that the expressions for 
the form factors obtained in PQChPT with N sea = 3 and 
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m a = mi) = m sea , indeed agree with the ones obtained 
in full ChPT with m u = m d = m s (for the full ChPT 
expressions see Ref . 0) • 



Form factors in the chiral limit 

In this section we make several important remarks con- 
cerning the chiral behavior of the form factors at some 
fixed (albeit small) v-p. We focus on the situation in 
which the light pseudoscalar meson consists of degener- 
ate valence quarks, i.e. My = 4/iora va i.- We distinguish 
the following three cases. 



(1) Expansion of (Sf p , v ) Loop (HHJ) for m sea and 
fixed nonzero m va i., results in a linear term in ra sea , 
but without logarithmic terms, i.e. 



C p ' v + CYM 2 S 



(19) 



where Mj 



4/io^sea- The coefficients Cq' v , Cf' v 



are functions of v-p and My, with ellipses repre- 
senting higher terms in the expansion. This re- 
sult is helpful for the future lattice simulations 
with 7V sea 7^ 0. Namely, it suggests that a lin- 
ear extrapolation of the form factors in m sea from 
the directly accessible sea quark masses down to 
m sea —> 0.5 x (m u + md) phys ' is not modified by the 
chiral logarithms, provided the valence quark mass 
is kept fixed. In the lattice studies, the constants 
Co and C2 are then to be obtained by fitting the 
data to Eq. JTJJ) . We will scrutinize this point more 
in the next section. 

(2) If, on the other hand, one keeps the sea quark mass 
m sea nonzero and studies the limit m va i. — > 0, the 
chiral logarithms appear. In addition to the ordi- 
nary chiral logarithms, i.e., of the form My ln(My), 
one also picks the quenched logarithmic divergences 
M|ln(M^). In this limit 



1 



^Loop 

h ~ (47r/)2AT S( 



(v-pY 
^M|ln(M^ 



1 + 3g 2 ) My In (My) 



^Txg^—^My 
v-p 



r Loop 



+ 3 5 2 



(20) 



(47T/) 2 iV s( 



c'l 



{2Ml - M 2 s )\n{Ml) 
-C1MI + ... , 



where the coefficients C%% are functions of v-p and 
Ms- We remind the reader that My = 4/io^vai.- 

(3) The case with m va i. = m sea is actually the case of 
full, unquenched, QCD. The form factors in the chi- 
ral limit, m va i. = m sea —> 0, are finite and coincide 
with the ones derived in the standard (unquenched) 



ChPT with N sea mass-degenerate quarks. In this 
limit the chiral behavior is 



J v,p 



Loop 



1 



(4tt/)2 



(l + 3g 2 ) (1 ^ J) M^{M>) 



c" v ' p + c'TmI 



• (21) 



Notice that the situation for the B — > K transition is 
qualitatively very similar to the B — > tt case, described 
above. 



IV. RESULTS ("PHENOMENOLOGY") 

As we mentioned in the introduction, the current 
quenched lattice simulations of the B — » tt matrix ele- 
ment are performed with the light quark masses m q > 
mP hys /2, where mf ys is the physical strange quark mass. 
This feature is likely to remain as such in the forth- 
coming (partially) unquenched lattice studies. Moreover, 
the first partially quenched lattice simulations will most 
probably be performed with 7V sea = 2. For that reason, 
in what follows, we consider the case of two degenerate 
sea quarks (A^ sea = 2). 

Once the partially quenched lattice QCD results be- 
come available, the NLO chiral expressions for B — > tt 
form factors (provided in Sec. lIII|) can be used to extrap- 
olate from the light quark masses directly accessed in the 
lattice simulations down to the physical w/d-quark mass. 

One way to proceed is to match the chiral expres- 
sions (fi7)|) with the lattice data at some intermediate val- 
ues of (m va i.)M and (m sea )M, at which the lattice data 
are used to fix the low energy constants (cf. discussion 
in Sec. 6.2 of Ref. 7). From that point down to the chiral 
limit, the extrapolation is made by using such determined 
constants, plus the coefficients of the chiral logarithms 
predicted in PQChPT. The matching procedure is needed 
to make contact of the pronounced linear light quark 
mass dependence of the (quenched) lattice data with the 
NLO expressions derived in ChPT, thus guiding the ex- 
trapolation to the physical pion mass. 4 It is, however, 
not clear at which point the above-mentioned matching 
should be made, i.e., that for the masses lighter than that 
used in the matching procedure one can trust the chiral 
perturbation theory. Is it mf ys , mP hys /2, or mP hys /4? 
Clearly, depending on what we choose for (m va i.)M and 
(^sea)M, from which we include the logarithmic terms 
in the extrapolation, we will get different results for the 
physically interesting form factors. Moreover, since the 
lattice data show linear dependence on quark masses this 



4 The linear dependence of the B —> tt form factors in the light 
quark mass is observed for the quenched data as well as for the 
B transition matrix element in both quenched and partially 
quenched studies. 
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also means that the variation of matching point will in- 
troduce a larger variation of the extrapolated values, if 
the ChPT dependence is very nonlinear 0, If in- 

stead the behavior of the form factor predicted by ChPT 
is close to linear, the precise point where we match ChPT 
expressions to the lattice data will not matter at all, as 
long as both the size and the slope (the leading and the 
NLO coupling of ChPT expression) are matched to the 
lattice data. It is thus far more reasonable to look for 
the strategy to extrapolate in m va i. and m sea , in which 
the ChPT expressions exhibit almost linear dependence. 
This is precisely where the discussion made in Sec. IIIII 
"Form factors in the chiral limit" becomes important. 
Recall that we found that the dependence of B — > ir 
form factors on m sea (with other variables fixed) is linear, 
whereas the one on m va i. does include nonlinear terms. 
However, one can suppress the most dangerous chirally 
divergent term ra sea lnra va i. if the m sea is close to the chi- 
ral limit. This suggests that a fairly linear behaviour can 
be obtained if one first extrapolates in m sea , and then in 
m va i.. It is this observation that we will elaborate more 
on in the present section. 

In the numerical evaluations we shall assume that the 
low energy constants appearing in PQChPT with N sea = 
2 are equal to their counterparts in the full ChPT with 
N = 3. In other words, we assume that the low energy 
constants depend weakly on the number of sea quarks, 
A/sea . In addition, and for an easier comparison, we will 
take the same values for the low energy constants as in 
Ref.S a = 0.56(4) GeV 3/2 , g = 0.5(1), / = 0.13 GeV, 
fio = 1.14(20) GeV, L 4 = -(0.5 ± 0.5) x 10" 3 , and 
L 5 = (0.8 ±0.5) x 10 -3 . The counterterms L4 and L 5 are 
evaluated at the renormalization scale \i = 1 GeV, which 
is also the choice for \i in the loop integrals. For more 
details about the choice of parameters and the complete 
list of references we refer the reader to Ref. 7. For the 
other counterterms ^1,2, ^1,2, there are no available ex- 
perimental data or any indications from the lattice data 
and we will take them to be zero. We have verified that 
the impact of the variation of these constants on our con- 
clusions is insignificant: (i) the loop corrections depend 
only on g and are not sensitive to variations of g in the 
range g = 0.5(1), (ii) the counterterms have negligible 
effect on the B — > ir form factors, while they can mod- 
ify the magnitude of B — > K form factors but not their 
dependence on v-p [Eq. JE}]. For the physical strange 
quark mass we take mf E (2 GeV) « 90 MeV 0. 



A. B — > 7r transition 

We work in the isospin limit and set the masses of the 
pion valence quarks to be equal, m a = = m va i.- To 
examine the chiral behavior of the form factors we will 



vp=0.2 GeV 




fval. 



= m va i./mf ys , 



phys 



(22) 



0.2 0.4 0.6 0.8 



sea 



FIG. 3: Comparison of complete f p ^ n expressions in the 
partially quenched N se & = 2 theory and the linear extrapola- 
tions. Here v-p = 0.2 GeV and r va i. = 0.2,0.5,1.0 are held 
fixed, while r sea is varied. Linear fit is made to the PQChPT 
expressions from Appendix B in the range Tsea g [0.5,1.0], 
and then extrapolated to the chiral limit. The vertical line 
marks the physical value r sea = 1/25. 



dimensionless parameters defined with respect to the 
physical strange quark mass. In the lattice studies, low- 
ering the sea quark mass, m sea (i.e., r sea ), is numerically 
costlier than lowering the valence quark mass, m va i. (i.e., 
r va i.). Fortunately, as discussed in the previous section, 
the dependence of f PlV (v m p) on r sea is expected to be lin- 
ear, so that even larger extrapolations in m sea may still 
lead to reasonably small extrapolation errors. This is il- 
lustrated in Fig. 03 for the case of the form factor f p: with 
v-p = 0.2 GeV. We see that for small r sea < r va i. and 
^se a —> 0, the form factor is indeed a linear function of 
r sea , as already stated in Eq. JTJJ). That behavior gets 
modified when r sea > r va i., resulting in the smooth forms 
that are very close to linear. If we linearly extrapolate the 
behavior of the form factors from the range of sea quark 
mass r sea G [0.5, 1] down to the chiral limit, we observe 
only a small discrepancy with respect to the complete 
expression for f p (of about 5 %). The same observation 
holds for the form factor f v , with the discrepancy be- 
tween the extrapolated and the result of the complete 
expression being approximately halved with respect to 
the one estimated for f p . As it can be observed from 
Fig. 03 the uncertainty induced by the extrapolation in 
^se a gets larger as the valence quark mass is decreased. 
However, for the quark masses that are realistic in the 
current lattice studies, r va i. > 0.5, the mentioned uncer- 
tainty is below the 5 % level. 

The other possibility, i.e., to extrapolate f PyV in r va i. 
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vp=0.2 GeV 




0.2 0.4 0.6 0.8 
val. 



FIG. 4: The situation similar to what is plotted in Fig. |3| 
but with r se a and r va i. exchanging the roles, i.e., r va i. is varied 
for three fixed values of r sea = 0.2, 0.5, 1.0. 



for a fixed value of r sea , is illustrated in Fig. 0] We see 
that in this case the difference between the linear ex- 
trapolation and the complete expression is much more 
pronounced, which reflects the presence of the quenched 
chiral logarithm ex r sea ln(r va i.), and hence the effect is 
larger for the heavier sea quark masses (cf. Eq. ((20)) ). 
Similar observations to the ones shown in Figs. [3| and 
are also valid for the form factor f v . 

An interesting observation is that the dependence of 
f PyV on ryai. (with r sea and v-p fixed) is such that the 
form factors at the physical point r va i. —1/25 are larger 
than what one gets from the linear extrapolation. On the 
other hand, when r va i. and v-p are fixed, the form factors 
at the physical point r sea ~ 1/25 are smaller than the 
linearly extrapolated ones. 

Therefore by first performing a linear extrapolation in 
r sea , and then in r va i. will not only allow one to avoid the 
spurious quenched logarithms, but it will also produce 
a mutual cancellation of the errors induced by the two 
chiral extrapolations. Indeed, the errors of the two con- 
secutive extrapolations, as shown on Fig. [51 are strikingly 
small. Performing linear fits to the N sea = 2 PQChPT 
expressions for r sea , r va i. £ [0.5, 1] leads to the errors be- 
low 5%, for both f p (v-p) and f v (v-p) form factors, and 
for a large range of v-p G [0.05,0.75] GeV. If, however, 
the chiral extrapolation is made linearly by keeping the 
valence and sea quark masses equal, r va i. = r sea , the re- 
sulting error is in the 20% range due to the explicit chiral 
logarithmic corrections given in Eq. (|2T)h 

For the reader's convenience, in Fig.[5J we also plot the 
result of the fully unquenched theory with N = 3. We 
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FIG. 5: Solid lines show f^^ 71 (v-p) in the partially quenched 
theory with A^ sea = 2 at physical point r sea = r va i. = 1/25. 
The dashed lines are obtained by first linearly extrapolating 
then in r va i. to the physical point, as explained in the 
text. For an easier comparison we also plot the same form 
factors obtained in ChPT with N = 3 by a dotted line (using 
the same values for low energy constants as in PQChPT with 

iVsea = 2). 



observe a pleasant feature that the shapes of the form 
factors in these two theories are very similar. Notice 
also that the results obtained in PQChPT with 7V sea = 
2 are systematically lower from the ones obtained with 
ChPT TV = 3 by 10% - 20%. It should be stressed, 
however, that this discrepancy might be artificial and 
merely a consequence of our assumption that the low 
energy constants do not change significantly when N = 
2^3. More information on the low energy constants in 
the theory with N = 2 is needed to get a clearer picture 
on the significance of the differences depicted in Fig. 

Another lesson comes from Fig. where we plot the 
rvai. dependence of the form factors by using the com- 
plete N sea = 2 expressions, derived in PQChPT, and by 
linearly extrapolating in r sea as discussed above. We see 
that the error due to the use of linear extrapolation of f p 
in r sea is fairly small at large r va i., but then increases as 
^vai. approaches its physical value. Reducing r va i. below 
^vai. < 0.5 in future lattice studies will therefore improve 
on the simple consecutive linear extrapolations (first in 
r sea , then in r va i., for r sea , 

fvai. £ [0,5,1.0]) only if, to- 
gether with smaller r va i., also the smaller values of r sea 
are reached in the simulation. For the high precision 
studies with both r sea and r va i. well below 0.5, it might 
even be more advantageous to use extrapolations with 
^sea = ^vai. (cf- Fig. left), contrary to what we advo- 
cated for more realistic studies (at present), namely with 
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FIG. 6: Solid lines show the r va i. dependence of f p ^ n 
[Eq. (|14f) ] in the partially quenched theory with A^ sea = 2 
and at r sea = 1/25, v-p = 0.2 GeV. The dashed lines show 
the form factors linearly extrapolated from r sea G [0.5, 1.0] 
down to physical point r sea = 1/25. The dotted lines refer to 
the dependence on r va i. while keeping 



^sea? ^val. ^ 0.5. 

Finally, we consider the soft pion theorem, which states 
that the ratio Fo/(/b//tt) —> 1 as v-p and quark masses 
go to zero in full ChPT. We verified that this is satisfied 
also in partially quenched theory in the limit in which 
v-p — > 0, and m sea = m va i. —> 0. Numerically, the ratio 
stays within 0.95 < [-Fo/(/b//tt)] < 1-25 in the ranges of 
v-p < 0.2 GeV and < r va i.,r sea < 1, but it diverges for 
m vaL -> if m sea / m vaL . 
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FIG. 7: Behavior of the SU(3) breaking ratio R PjV = 
fp,^T K I fp,^ 7 " physical values of quark masses) as func- 
tions of v • p. The prediction of PQChPT with A^ sea = 2 is 
denoted by the solid lines. The dashed curves are obtained 
after consecutive extrapolations, first m Tsea and then in r va i. 
(as explained in the text), while the dot-dashed curve is the 
prediction of the full ChPT with N = 3. 



B. Ratios f*? K /f*f 

Instead of repeating the previous discussion for the 
B — > K form factors, we investigate the ratios R PiV — 
fp,v* K /fp,v*^- Our kaon is composed of valence quarks 
with one mass fixed at the physical strange quark 
mass and the other varied, i.e., Mj^ = 2/j J Qm^ hys (l + 
r*vai.)> whereas the pion valence content reads, M% = 
4/io^s hys r va i.- We also ignore the difference between the 
poles in fp^ K and in fp^ n , i.e., we take tub* = tub*- 

We again observe the linear dependence of the ratio 
R PjV under varying the light sea quark mass, m sea , while 
keeping fixed values of m va i. an d v-p. As in the previous 
section we see that the variation of m va i. m the chiral 
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limit at fixed m sea exhibits both standard and quenched 
chiral logarithmic dependences. Therefore, like in the 
B — > 7r case, a consecutive linear extrapolation in r sea and 
then in r va i. leads to very small errors of extrapolation, 
namely below 5% (see Fig. 0. 

In Fig. [7| we observe that the dependence of R p on v-p, 
as inferred from the PQChPT with N sea = 2, differs from 
the one obtained in the full ChPT with TV = 3 at the 10% 
level. Note that this feature does not change under the 
variation of C(p 3 ), 0(p 4 ) low energy constants given in 
Eqs. Q and (0. In addition, the v-p dependence of R p in 
partially quenched theory does not change significantly 
when varying r sea in the range r sea G (0, 1). It is impor- 
tant to stress that the apparent disagreement between 
^PQChPT anc j ^>ChPT ag f unc tions of v-p disappears after 
replacing N sea = 2 by N sea = 3 in our PQChPT expres- 
sions. Finally, it should be mentioned that the absolute 
values of R PyV (but not their v-p behavior) may be fur- 
ther affected by the uncertainties in H\ and L5, because 
these constants enter the expressions for R PiV , multiplied 
by rn 2 K . 



a result, the above procedure leads to results which are 
different by no more than 5% with respect to the com- 
plete prediction of PQChPT with N sea — 2. This makes 
a strong call for the partially quenched lattice compu- 
tations of the heavy — > light transition matrix elements 
with 7V sea = 2. 

Regarding the comparison of the form factors obtained 
in PQChPT with A/g ea = 2 and the ones obtained in 
ChPT with N = 3, we notice that the results obtained 
with N = 2 are systematically lower than their N = 3 
counterparts. The difference is in the ball park of 10% — 
20%. At this point it is not clear whether this difference 
is a realistic estimate of the present approach, or merely 
a consequence of our lack of knowledge of the values of 
the low energy constant in the theory with N = 2. The 
above estimate of 10% — 20% is obtained by assuming the 
low energy constants to be independent of the number of 
flavors. 
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V. CONCLUSIONS 

In this paper we presented the expressions for the 
heavy — > light pseudoscalar meson form factors as de- 
rived at NLO in PQChPT and in the static limit of 
HQET. This calculation is useful in the perspective of the 
QCD simulations on the lattice, since the most wanted 
B — > 7r transition form factors are expected to be com- 
puted quite soon in the partially quenched QCD with 
7V sea = 2. The approach adopted in this paper is valid 
for small recoils (i.e., close to q 2 — > ^ ax ), the part of the 
phase space that is also accessible by the current lattice 
simulations. The main benefit of the expressions derived 
in PQChPT lies in the fact that one can disentangle the 
chiral behavior in the sea quark mass, m sea , from the 
chiral dependence in the valence light quark mass, m va i.. 
Our results clearly suggest a linear behavior of the form 
factors under the variation of the sea quark mass value 
m sea , when the valence quark mass m va i. is kept fixed. 
On the other hand, if the sea quark mass m sea is fixed to 
some value that is directly accessible from the current lat- 
tice studies with N sea = 2 (e.g. r sea > 0.5), the form fac- 
tors exhibit both the standard and quenched chiral loga- 
rithms when approaching the physical point, r va i. ~ 1/25. 
For that reason, we propose a strategy of computing the 
B — > tt(K) semileptonic form factors on the lattice (in 
partially quenched theory with 7V sea = 2 degenerate light 
quarks) that requires a computation to be performed at 
several values of the sea quark mass with the valence 
quark mass held fixed. The extrapolation of F +5 o(<? 2 ) in 
the sea quark mass to the physical u/d quarks can then be 
made linearly at fixed valence quark mass. This should 
then be followed by an extrapolation in the valence quark 
mass. In this way the quenched logarithm is avoided since 
the pathological term is of the form m sea ln(ra va i.)- As 
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APPENDIX A: CHIRAL LOOP INTEGRALS 

In this appendix we list the dimensionally regularized 
integrals encountered in the course of calculation. For 
more details, see Ref. 20| and references therein. 

tft J ( 27 r)4-V-m2 = 16^ Jl(w) ' 

1(1 J (2tt) 4 - (q 2 - m 2 )(q-v - A) " 16tt 2 A j ' 

(Al) 

where 

1 1 (rn) = rn 2 ln^^-^j — m 2 A , 

I 2 {m, A) = -2A 2 ln(^) - 4A 2 f(^) + 2A 2 (1 + A) , 

(A2) 

where A = 2/e - 7 + ln(47r) + 1. The function F(x) 
has been calculated in Ref. |2j, for both the negative and 
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positive values of the argument: 



vT 



7T 

tan 

2 



\x\ < ^ 



— In [x + \/^ 2 — 1^) 



The functions Ji(m, A), J 2 (m, A) differ from the ones 
in Ref. [22 by the last two terms in Eq. (|A6|) which are 
,of 0(m 2 ,A 2 ). These additional (finite) terms originate 
rom the fact that r\^ v is the (4 — e)-dimensional metric 
tensor. 



N >i. 

(A3) 



In addition to the integrals (|A1|) . one also needs the fol- 
lowing two: 



16tt 2 



(27r) 4 - e (<? 2 - m 2 )(q-v - A) 
[J 2 (m,A) + Ji(m)] , 



(A4) 



s 

1 

16tt 2 



d 4 " 



(27r) 4 - e (<? 2 - m 2 )(q-v - A) 



A[Ji(m,A)<^ + J 2 (m,A)^^] , (A5) 



with 



Mm, A) =(-m 2 + ?A 2 ) In fej + | (A 2 



^ A 2/i i A\ i ^_2/ i oA\ i 



- -A^(l + A) + -m\2 + 3A) + -m 



J 2 (m, A) =(2m 2 - ^A 2 )ln _ !( 4 A 2 
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-A 2 (l + A)--m 2 (l + 3A)--m 



9 

(A6a) 
(A6b) 



APPENDIX B: LOOP CONTRIBUTIONS TO 
FORM FACTORS IN PQCHPT 



In this appendix we list expressions for the one-loop 
chiral corrections to the form factors f PjV [Eq. ([T£])] 
in partially quenched ChPT (the form factors for fully 
quenched ChPT and full ChPT are collected in Appendix 
B of Ref. 7). The formulas apply to the B a — > 
transition with B a ~ bq a and P\y a ~ qbQa- They are 
expressed in terms of the pseudoscalar meson masses 
M I,b,s = 4 Mo^a,6,sea, M 2 b = 2/i (ra a + mb), and 
M 2 ShS = 2/io(^a,6 + ^sea)- We list the nonzero contri- 
butions to /pfj, with the superscript "J" corresponding 
to a label of the diagram shown in Fig. 



J 



*4 13) = 

and 

Sfi 4) 



(14) 



(4tt/) 2 I 



Ji(M bS ,vp) - — ^M 3 S 



TV. 



1 f l + (Mi-Ml) ° 



1 27T 

Ji(M b ,vp) -M 3 

3 



i 



i 



(4tt/)2 AT sea {Ml - Mi 



(Ml Ml) [ ,h(M a ,v.p) - M a 3 ) - {Ml - Ml) (.h(M b ,v.p) - ±-*± M 3 



6(4tt 
1 



■^2 |iV S ea(/i(M oS ) + /i(Mi,s)) - ^- (Gi(M a , M s ) + Gi(M 6 , M s )) + -j^—H\{M a , M b , M s ) J , 



2(4tt/) s 



2(4tt/) 2 



N sea h{M bS ) 



-G!{M b ,M s ) 



Ii{M bS ) + 2I 2 (M bS ,vp) 



1 



TV sea 
2 



/Vs, 



6(4tt. 



t/) 2 { 



Ji(M oS )+Ji(M 6S ) 



iV s< 



H 1 (M a ,M b ,M s )-G 1 (M b ,M s ) 
H 2 (M a , M b , M s , v-p) - G 2 {M b , M s , v-p) 

1 



G 1 {M a ,Ms) + G 1 (M b ,M s ) 



ffi(M a ,M 6 ,M s ) 



(B2) 
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The one loop chiral corrections to the wave function renormalization factors Zb p are 



£ Z Loop = 

P ab 3 ( 47r 



r\o \ sea 

-NsMMas) 



h{M aS ) + h{M bS ) 
1 



N„ 



Gi(M ,M s )| . 



ffi(M a ,M b ,M s ) 



1 



G 1 (M a ,M s ) + Gi(M 6 ,M s ) 



(B3) 



The functions Ii(ra), him), and Ji(m, A) are given in 
Appendix A. In addition, we have used the abbreviations 

Gi (M, M S ) = ^ [(M 2 - Mj)/! (M)] , (B4) 

Hi(M OJ M 6 , M 5 ) = M2 ^ M2 [(M a 2 - M|) Jx (M a ) 

- (M^-Ml)/^)] , (B5) 

and similarly for G 2 (M, M s , v-p) and H 2 (M a , M b ,M s , v- 
p), that are obtained by simply replacing ii(M) — > 
l2(M,v-p) in the above expressions. Note that in the 
limit of degenerate valence quarks, M a — > M5, we have 

, lim ffi(M a , M b , Ms) = Gi(M a , M s ) , (B6) 
lim H 2 (M a ,M b ,Ms,v-p)=G2(M a ,M s ,v-p). (B7) 



We reiterate that the mass differences (A) among 
B*, £? s , and B J mesons have been consistently neglected 
in the loops, since v-p > A*, for the B — > 7r, if tran- 
sitions. This induces a spurious singularity in the ex- 
pression for the diagrams (7), at v-p — > 0. To handle 
such singularities we followed the proposal by Falk and 
Grinstein 23] and resum the corresponding diagrams and 
then simply subtract the term that would renormalize the 
5*-meson mass. 
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